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We consider general metric f{R) theories of gravity by solving the field equations in the presence of 
a spherical static mass distribution by analytical perturbative means. Expanding the field equations 
systematically in 0{G), we solve the resulting set of equations and show that f{R) theories which 
attempt to solve the dark energy problem very generally lead to ^ppn ~ 1/2 in the solar system. 
This excludes a large class of theories as possible explanations of dark energy. We also present the 
first order correction to "/ppn and show that it cannot have a significant efi'ect. 



I. INTRODUCTION 

The dark energy problem remains central in modern 
day cosmology. Since the matter only, homogeneous uni- 
verse within the framework of general relativity is in con- 
flict with cosmological observations, the assumptions be- 
hind this model have been questioned. The most popular 
modification is to consider a universe filled with other, 
more exotic forms of matter, the cosmological constant 
being the leading natural candidate. Other ways to tackle 
the dark energy problem are then to relax the assumption 
of homogeneity or modify the theory of gravity. 

In recent years, a particular modification of gravity, the 
f{R) gravity models that replace the Einstein-Hilbert ac- 
tion of general relativity (GR) with an arbitrary function 
of the curvature scalar (see e.g. P, S i, II i, 11 0, 11 
and references therein) have been extensively studied. 
Naive modification of the gravitational action is not with- 
out challenges, however, and obstacles including cosmo- 
logical constraints (see fssl . Isol . l40j and references 
therein), instabilitie s ll]] . [l3j . solar system constraints 
(see e.g. [T3.[l5l[l^l2l| and references therein) and evo- 
lution large scale perturbations [13, [H, need to be 
overcome. In addition, a number of consistency require- 
ments need to be satisfied (see e.g. and references 
therein) . 

One of the most direct and strictest constraints on any 
modification of gravity comes from observations of out 
nearby space-time i.e. the solar system. This is often 
done by conformally transforming the theory to a scalar- 
tensor theory and then considering the Parameterized 
Post- Newtonian (PPN) hmit [23.,^] (see also [25, 26] for 
a discussion). The question of validity of the solar system 
constraints f{R) theories has been extensively discussed 
in the literature and not completely without controversy. 
The opinions on the viability of f{R) theories have been 
divided from more or less skeptical [13, HI, HI, 1131 to 
approving [3l|, depending on the point of view of the 
author. 

The essence of the discussion has been the question 
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of validity of the Schwarzschild-de Sitter (SdS) metric as 
the correct metric in the solar system. The SdS metric 
is a vacuum solution to a large class of f{R) theories of 
gravity but due to the higher-derivative nature of metric 
f{R) theories, it is not unique. Other solutions can also 
be constructed in empty space, in the pre sence of matter 
and in a cosmological setting (see eg. [lOl. IbsI. Ib^I). 

In hght of recent literature [13,[li,[l^, the validity of 
the solar system constraints has become clear and it is 
now understood that the equivalent scalar-tensor theory 
results are valid in a particular limit that corresponds to 
the limit of light effective scalar in the In terms of the 
/(i?) theory, this is equivalent to requiring that one can 
approximate the trace of the field equations by Laplace's 
equation [2^]. As a result, the often considered R — fi^/R 
theory [J (the CDTT model) is not consistent with the 
Solar System constraints in this limit, if the 1/R term is 
to drive late time cosmological acceleration. 

In [l^l the CDTT model was considered by lineariz- 
ing around a static de Sitter spacetime and solving the 
trace equation in terms of R(r), resulting in a spacetime 
outside the star where 7 = 1/2. This result was then 
generalized for a general f{R) theory in [2^ by study- 
ing the space-time outside a spherical mass distribution 
and expanding f{R) in terms of a perturbation in R. 
Again solving the trace equation leads to an outside so- 
lution with 7 = 1/2 as long as the effective scalar mass 
is light. A somewhat different approach was followed in 
1^], where the trace equation was first written in terms 
F{r) = df /dR in the perturbative expansion. Solving the 
trace equation then leads to 7 = 1/2 outside the star. 

In this paper we follow the latter approach by view- 
ing F along with the metric as independent functions. 
By expanding all quantities in G and solving the result- 
ing equations inside and outside the star for a general 
f{R) theory, we find that generally, jppN = 1/2 + 0{G) 
outside the star and the scalar curvature is ©(G^) ev- 
erywhere. We also idenfity the first order correction to 
IPPN and show that it cannot have a significant effect. 
Only if initial conditions inside the star arc fine-tuned 
such that the scalar curvature follows the matter density 
like in GR [U, [s^ can these bounds be evaded. 
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II. f{R) GRAVITY FORMALISM 

The action for f{R) gravity is (c = 1) 

5 = 1 d^x^[j^f{R) + . (1) 

The field equations resuhing in the so-caUed metric ap- 
proach arc reached by variating with respect to g^^: 

F{R)R^,-]^f{R)g^,-V^y,F{R)+g^,nF{R) = SnGT^.^ 

(2) 

where is the standard minimaUy coupled stress- 
energy tensor and F{R) = df /dR. 



Contracting the field equations and assuming that we 
can describe the stress-energy tensor with a perfect fiuid, 
we get 

F{R)R - 2/(i?) + 3nF(i?) = SnGip - Sp). (3) 

In this letter we consider spherically symmetric static 
fiuid configurations and adopt a metric, which reads in 
spherically symmetric coordinates as 

ds^ = B{r)dt^ - A{r)dr^ - r^d0^ - sin^ 0d(p'^ . (4) 

By taking suitable linear combinations of the field equa- 
tions they can be written in the following form: 
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where prime indicates a derivation with respect to r, ' = 
d/dr and we have written / and F as functions of the ra- 
dial coordinate r expect in combination 2f{R) — RF{R), 
which we will expand in terms of curvature R. 
The corresponding equation of continuity is 

P'jr) ^ IB'jr) 
p{r)+p{r) 2B{r)- ^' 

When pressure is negligible, it is easy to see that B must 
be a constant. This is, however, not acceptable and there- 
fore an adequate perturbation expansion is needed. 

III. PERTURBATIVE EXPANSION AND ITS 
SOLUTIONS 

We expand the metric as well as F with G as an ex- 
pansion parameter: 

A{r) = l + GAi(r)+0(G2), 

B{r) = Bo + GBi{r)+0{G^), 

F{r) = Fo + GFi{r)+0{G^), 

p{r) - po + Gpi(r)+0(G2). 

Note, that we consider the density profile p{r) to be a 
fixed function and also that Bq and Fq are constants. 
From the expansion of A and B, one can also read out 
an expansion for R: 

R^ Ro + GRi+0{G^). (9) 



From the equation of continuity we see that at 0{G^) 
pressure is constant and exactly zero, po = 0, simply 
because it vanishes in empty space. Therefore pressure 
effects are always ©(G^) and do not contribute to the 
0{G^) expansion. 

The 2/ — FR term in the third field equations is cru- 
cial in determining the behaviour of the solution. In gen- 
eral, for a f{R) dark energy model, this term is neg- 
ligible and can be omitted, at least in the first order 
approximation. This is demonstrated explicitly for the 
CDTT in model and discussed more generally in [35| . 
where it is argued that the non-linear term is completely 
negligible, barring fine tuning. This argument is easily 
understandable in a general model since in the vacuum 
2/ — FR ~ GpDE <C Gp for any stellar matter configu- 
ration. Note that this will in general be true also outside 
a stellar configuration as the dark matter will completely 
dominate over the cosmological term. Hence, in the trace 
equation, the non-linear terms can be dropped, unless 
the initial conditions are fine-tuned. We will return to 
the fine-tuned condition, or the Palatini limit [2ll, [ssj . 
later. 

More formally, the same conclusion can be confirmed 
by using an expansion in G for the non linear-terms as 
well: 2/(i?) - = 2/(i?o) - F(i?o)i?o + {F{Ro) - 

F'{Rq)Rq)Ri+0{G'^) Evidently, the expansion point Rq 
has to be such, that it corresponds to the correct back- 
ground of the theory, i.e. 2/(i?o) - F{Ro)Ro = 0. Then 
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expanding up to first order in G, the field equations are 
FoA[ , FoB[ 



r Bor 
FoA, FoA[ , FoB[ F[ , F^B'{ 

r 



2r 



2B, 



F[' ^ Snp, 

Snp, (10) 







2Bor 

hRi + ^+3F[' = -8np, 
r 



where Ii = F(Rq) — F'{Rq)Rq is a constant and Fq = 
FiRo). 

The set of equations (fTO]) can be straightforwardly 
solved leading to 0{G) functions: 
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Inserting this solution back to expression of the curvature 
scalar, we find that i?i = 0, i.e., R is 0{G^). It is crucial 
that in deriving the solution (jlip. we have assumed that 
A, B, F are regular at the origin. 

The PPN-parameter is now straightforwardly calcula- 
ble: 
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It is easy to see that, at the boundary of the star 
IPPN ^ 1/2 + 0{G). This behaviour was also observed 
in numerical studies [s^ [s^ . From the first order cor- 
rection one can furthermore conclude that if one wishes 
corrections to be effective at zeroth order, Fq needs to be 
of order G. However, looking at the continuity equation, 
Eq. dSD, we find that 



ry ^~pm{r)+0{G^). 

6 Po 



Comparing this with the Newtonian result. 



(16) 
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AGpm{r), we see that if Fq ~ 0{G), the effective New- 
ton's constant is orders of magnitude larger than the one 
in Newton's theory (or GR), resulting in stars with a 
completely different mass to radius relationship than the 
one observed. Furthermore, from the continuity equa- 
tion, we can read that unless Fq « 4/3 to a high pre- 
cision, a star with the same density profile, and hence 
total mass, will have a different radius than in GR. This 
behaviour was already observed in [3^ . 

In general the results described in this section will 
apply even when 2/ — FR ^ R, i.e. when f{R) = 



R+ciR^. Because in the approximation described above, 
R ^ 0{G^), it is easy to see that that this term will play 
no role in the trace equation. Similarly for higher order 
terms in R. One can avoid the constraint only if F has 
no G order correction. In this case, the DF, term is neg- 
ligible in trace equation and we recover the GR results, 
or the Palatini limit [2l|, [3^ . Alternatively, if one relaxes 
the regularity constraint of the metric at the origin, one 
can also avoid the constraint as demonstrated in [36| for 
the CDTT model. 



A. Recovering the general relativity 

In the Palatini limit, where the trace-equation is sim- 
ilar than in the Palatini formalism, the theory is fine- 
tuned so that 2/ — FR fa R fa —SirGp throughout (see 
^35*1 for a numerical example). This is the mechanism 
that allows one to construct solutions that are consistent 
with solar system observations In the Palatini limit, 
the field equations read as 
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The "fppN parameter is easily calculable: 



jPPN ~ 1 - r— + 0{G). 
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Therefore, in this limit ^ppN 1 at the surface of the 
star. However, as shown in 35] for the CDTT model, 
this limit can be unstable time leading to the Dolgov- 
Kawasaki instability 11]. Stable theories are considered 
in j21| and studied analytically in [svj . 



IV. DISCUSSION AND CONCLUSIONS 

In this letter we have considered a general metric 
f{R) theory in the presence of matter by analyzing the 
field equations by perturbative means in linear order in 
the Newton's constant G. We have shown explicitly 
that for a typical star, any modification of gravity from 
GR will naturally lead to physically unacceptable value 
IPPN — 1/2. This places a very strong constraint on any 
,f{R) theory, in particular when acting as a dark energy 
candidate. Furthermore, even if the gravity theory is not 
motivated by cosmology, but by other arguments, such 
as quantum gravity, the presence of non-linear terms can 
still lead to a space-time inconsistent with observations. 

In this order of perturbation theory we can recover 
the observationally acceptable space-time, only when 
F — df / dR has no order G correction. Such a constraint 
indicates fine-tuning in the initial values of the solution 
so that one remains in the high curvature limit, R ^ Gp 



4 



throughout. However, the stabihty of such a fine-tuned 
solution may be p roblematic [35| . although possible to 
obtain [Hi ^ [H] . 

Since our analysis is of order 0{G^), further study on 
the system, in particular second order perturbations in 
G, may affect the conclusions. Indeed, our analysis shows 
that the first order perturbation theory is essentially in- 
dependent on the details of the underlying f{R) theory. 
The only piece of information used was the knowledge 
that there are higher order derivatives in the equations 
of motion, i.e. that the theory is not GR. New effects 
may appear in higher order perturbation theory, where 
finally the dependence on the functional form of f{R) 



should become evident. However, our results suggest that 
unless the solution is fine-tuned so that R ~ Gp through- 
out the mass distribution, a naive modification where a 
small correction is added to the Einstein-Hilbert action 
to solve the dark energy problem is not likely to pass the 
solar system constraints. 
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